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Abstract. 

Braunstein and Caves (1994) proposed to use Helstrom's quantum information 
number to define, meaningfully, a mètric on the set of all possible states of a given 
quantum system. They showed that the quantum information is nothing else than 
the maximal Fisher information in a measurement of the quantum system, maximized 
over all possible measurements. Combining this fact with classical statistical results, 
they argued that the quantum information determines the asymptotically optimal rate 
at which neighbouring states on some smooth curve can be distinguished, based on 
arbitrary measurements on n identical copies of the given quantum system. 

We show that the measurement which maximizes the Fisher information typically 
depends on the true, unknown, state of the quantum system. We close the resulting 
loophole in the argument by showing that one can still achieve the same, optimal, rate 
of distinguishability, by a two stage adaptive measurement procedure. 

When we consider states lying not on a smooth curve, but on a manifold of higher 
dimension, the situation becomes much more complex. We show that the notion of 
"distinguishability of close-by states" depends strongly on the measurement resources 
one allows oneself, and on a further specification of the task at hand. The quantum 
information matrix no longer seems to play a central role. 



PACS numbers: PACS numbers: 03.65.Bz, 03.67.-a 
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1. Introduction 

Braunstein and Caves (1994) have clarified the relation between the classical Fisher 
expected Ínformation number i(9), for the unknown parameter 6 of a probability 
distribution p(x;9), and the analogous concept of expected quantum Ínformation 1(9) 
for a quantum system in state p = p(9) on some Hilbert space. They showed that 
1(9) is the maximal Fisher Ínformation i(9; M) in the distribution of the outcome of 
a measurement M, over all measurements of the state. Thereby they supplied a new 
proof of Helstrom's (1967) quantum Cramér-Rao bound: no unbiased estimator of 9, 
based on any measurement, has variance smaller than I(9)~ 1 . Recali that the classical 
bound states that no unbiased estimator of 9 based on the outcome of the measurement 
M has variance smaller than i(9; M) _1 . 

For n identical copies of a quantum system, and for n independent and 
identically distributed observations from a probability distribution, quantum and Fisher 
Ínformation are both n times the corresponding quantities for n — 1. By classical 
statistical theory, the quantum bound is therefore asymptotically achieved, as n — > oo, 
by the maximum likelihood estimator of 9 based on the outcomes of the measurement 
maximizing the Fisher Ínformation for n — 1, applied to each of n copies of the quantum 
system separately. 

In the present paper, we analyse the conditions for equality of the quantum and 
Fisher Ínformation. We show that in general there does not exist a measurement M 
such that i(9; M) = 1(9) for all 9 simultaneously, studying the pure state, spin-half 
case in detail. In that case the model describes a curve on the surface of the unit 
sphere, specifying the direction of the spin as a function of 9. We show that one has 
uniform attainability if and only if the curve is a segment of a great circle. We show 
how (in general) adaptive measurements still allow one to asymptotically achieve the 
quantum Ínformation bound for a scalar parameter, though not in the vector case, where 
the picture is rather complicated and the quantum Ínformation matrix inadequate to 
describe what is possible. 

In Section 2 of the paper we recapitulate some of the theory of classical and quantum 
Ínformation. Next, in Section 3, we specialize the conditions for attainability of the 
Ínformation bound, first to pure states, then further to spin-half models. Unless the 
model specifies a great circle, no measurement achieves the bound uniformly in the 
parameter 9. In Section 4 we explore the consequences of this result. We show that 
one can in effect achieve i(9; M n ) pa nl(9) for all 9 simultaneously, when we measure n 
identical copies of the quantum system in one joint measurement M n . This result gives 
support to Braunstein and Cave's interpretation of the quantum Ínformation number 
1(9) as a measure of statistical distinguishability between neighbouring quantum states. 
Finally we turn to the case when the parameter is a vector. Both quantum and classical 
Ínformation numbers have matrix generalizations, and inequality between them still 
holds, in the sense of positive semi-definite matrices. The inequality is sharp but however 
no longer attainable. For a completely unknown spin-half pure state we show that the 
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optimal rate at which one can distinguish between different states does not follow from 
the quantum information matrix in the way one would expect from analogy with classical 
Fisher information. Moreover it depends on some weighting of the different aspects of 
the states which one wants to distinguish. Major open problems remain, and the role 
of the quantum information does not appear to be primary. 

A preliminary version of this paper appeared as Barndorff-Nielsen and Gili (1998). 

2. Expected classical and quantum information 

On a given Hilbert space, consider a quantum state (density operator) p = p{6), which 
depends on an unknown scalar parameter 9. Consider also a generalised measurement 
(operator-valued probability measure, POVM) M with outcomes in a measurable space 
(X,A). Thus the outcome of a measurement of M on p is a random variable X taking 
vàlues in X, such that for each measurable subset A of X, i.e., for each A e A, we 
have Pr e {X G A} = tr &ce p(9)M(A). Suppose that M is dominated by a sigma-finite 
measure p on (X,A), i.e., for each A E A, 

M(A) = / m(x)p(dx) 
Ja 

where the operator m(x) is, for each x, nonnegative and selfadjoint, and f x m(x)p(dx) = 
1. (This is no restriction for finite dimensional quantum systems, for which one can 
always take p to be the measure defined by p(A) = traceM(A).) Under the domination 
assumption, the outcome X of a measurement of M on p has probability density, with 
respect to p, given by 

p(x;9) = trace p(9)m(x) . 

Under sufficent smoothness, the expected Fisher information number, for 9, from this 
measurement, is defined by 

i{9;M) = E e Í(9) 2 = I (Í(x;9)) 2 p{x;9)p(dx) 

Jx 

where 

1(9) = l(X;9) = \ogp(X;9) 
is the log likelihood and 

is the score function for 9. 

Now, let X = X(9) denote the symmetric logarithmic derivative of p with respect to 
9, that is, the self-adjoint operator given implicitly by 

p = |(pA + Ap). (1) 
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We call À the quantum score for 9. From the relation trace p = 1 one finds, by 
differentiating, trace pX = 0. The expected quantum ínformation number for 9 is defined 
by 

1(9) = trace pA 2 . 

Note triat this quantity is defined without reference to any particular measurement M. 

For future reference, we mention that when 9 is a vector parameter, the Fisher 
Ínformation matrix is defined in the obvious way, while the quantum Ínformation matrix 
has r/'-component ^trace p(AjAj + XjXi) where Xi(9) is the quantum score for 9i keeping 
the other components of 9 fixed. For completeness, we mention that there exist other 
generalizations of quantum Ínformation; see Yuen and Lax (1973), Belavkin (1976) and 
the books Helstrom (1976), Holevo (1980). 

Define 

Xq = X (M,9) = {x:p(x;9)=0} (2) 

and let X + be its complement. One can express the Fisher Ínformation i(9; M) in terms 
of the quantum score for p: 

i(9;M) = I p(x;9)- 1 (^tiace(pXm(x))) 2 p(àx). 
Jx+ 

This follows on noting that 

1(9) = p(x; 9) _1 trace pm(x) 

= p(x; 6>) _1 |trace((pA + Xp)m(x)) 
= p(x]9y 1 ^R,tia,ce(pXm(x)). 

The usual proof of the quantum Cramér-Rao inequality, see Helstrom (1976) or, 
for a more abstract and precise version, Holevo (1982), follows closely the lines of the 
usual proof of the classical bound: write down the unbiasedness relation, differentiate 
under the integral sign, and apply the Cauchy-Schwarz inequality; Holevo (1982) gave 
a more rigorous proof on the same lines. Braunstein and Caves (1994) noted how the 
quantum bound could be obtained from the classical bound together with their new 
inequality i(9;M) < 1(9) for all measurements M. The derivation is a chain of three 
inequalities and therefore leads to a set of three necessary and sufncient conditions for 
equality (though they did not notice the third). Before presenting the derivation we 
list the three ingredients. For a given x G X + , let A = m(x) 1 ^ 2 p 1 ^ 2 , B = m(x) 1 ^ 2 Xp 1 ^ 2 , 
and z = tr&ce(A*B). The the first inequality step uses the trivial ($l(z)) 2 < \z\ 2 
with equality if and only if $>(z) — 0. The second uses the Cauchy-Schwarz inequality 
|trace(A*_B)| 2 < tr&ce(A* A)tr&ce(B* B) with equality if and only if tr&ce(A*A)B = 
trace(A*B)A. The third inequality step, trace(M(A' + )ApA) < trace(pA 2 ), follows from 
the fact that M(X + ) = 1 - M(X ) where M(X ) > 0. The three ingredients are put 
together as follows: 

i(9;M) = I p(x;9)- 1 (^trace(pXm(x))) 2 p(dx) 
J x+ 
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trace( (m(x)2p2)* (m(x)2Ap2) ) (trace(pm(x))) 1 /j,(dx) 



Ix+ 

< / trace(m(x)ApA)/z(da;) 
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= tmce(M(X + )\p\) 

< trace(pA 2 ) = 1(9) . (3) 

With A, B and z = trace(A*£>) as above (depending on x), necessary and suficient 
conditions for equality at the first two inequality steps in (3) together are equivalent to: 
for p(dx) almost all x in X + , tiace(A* B) is real and A ocjk B, by which we mean A = rB 
or B = rA for some real number r. But if A oc R B then automatically trace(A*_B) is 
real. Thus we have equality in (3) if and only if the following two conditions are satisfied: 
firstly for p,(dx) almost all x in X + 

11 11 
m(x)z\p2 oc M m(x) 2p2 (4) 

and secondly, 

trace(M(A' )ApA) = 0. (5) 

Obviously a sufficient condition for (5) is that M(X ) = 0, and a sumcient condition for 
that is p(x; 9) > for all /i-almost all x. Braunstein and Caves remark that a sumcient 
condition for (4) is that each m(x) is proportional to a projector onto an eigenspace 
of A. In particular, if the measurement M is a simple (von Neumann) measurement 
of the observable A then (4) is satisfied. However this is not a necessary condition for 
attainability. Thus the obvious fact that, in general, A(#) varies with 9, does not show 
that there are no measurements attaining the bound (3) for all 9 simultaneously. We 
will do this by a further study of condition (4) in a special case. 



3. Attainability of the quantum information bound 

In this section we concentrate on models for pure states, p = \ip) (ip\ where = \ip{9)). 
In this case, the quantum score can be computed explicitly and the condition (4) 
simplifies. Define the (unnormalized) state \a) = Since p 2 = p, we have 

p = pp + pp. The defining equation (1) for the quantum score therefore telis us that 
A = 2p = \a) (ip\ + \ip) (a\. Now let |1) = \ip) and let |2) be a normalized orthogonal 
state such that \a) is in the subspace spanned by |1) and |2); write \a) = a\ |1) + a<i |2) 
where a± = (1 | a) and a 2 = (2 | a). (Note that all these definitions are relative to a 
given value of the parameter 9.) We find that 

pï\ = p\ = 2(3?ai)|l)(l] + a 2 |l)(2| 
and condition (4) reduces to 

m(x)5(2(ftai)|l)<l| + a 2 |2)(l|) oc R m{x)\ |1) (1| . 
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This can be again simplified, resulting in the condition 

i i 
a 2 m(x) 2 |2) (1| oc R m(x)2|l)(l| 

or equivalently 

i i 

a 2 m(x) 2 |2) ocr m(x)2|l). (6) 

For spin-half models, thus a Hilbert space of dimension 2, a further simplification 
occurs. One can take |2) = \ip) ± , forming an orthonormal basis (depending on 9) with 
|1) = \ip). /,From (6) it follows that if m(x) satisfies (4), it must have less than full rank, 
and hence in the two-dimensional case both it and its square root must be proportional 
(with real constants of proportionality) to |£) (£| for some normalized state |£) = 
A minor rewriting yields that (4) is equivalent to the statement: for p(x; 9)p(àx) almost 
all x, m(x) is proportional to a one-dimensional projector \£(x)) {£(x)\ satisfying 

<£|2)<2|a) oc K (7) 

We show that this algebraic condition has a simple geomètric interpretation. 
First note that from the definition of \a) and the fact that (ip \ ip) = 1 for all 9, 
it follows that 23? (a | 1) = (a | 1) + (1 | a) — 0, hence (a | 1) is purely imaginary. 
By multiplying |2) by a suitable phase factor, one can arrange that (a | 2) is real 
and (7) becomes (£ | 2) ocr (£ | 1). Note that 2p = |1) (a\ + \a) (1|. It follows that 
(1 | f, | 1) = = (2 | p | 2), while (1 | p | 2) is real. Hence 2p = r(|l) (2| + |2) (1|) for some 
real number r. Let a x = |1) (2| + |2) (1|, a y = -i |1) (2| +i\2) (1|, <r 2 = |1) (2| + |2) (1| 
be the Pauli spin matrices with respect to the basis |1) , |2). In this basis p = |(1 + a z ) 
and 2p = ra x . We can write |£) (^| = |(1 + a • a) where a = (a x , a y , a z ) is a unit vector 
in IR 3 . Note that (2 | ^) | 1) = a x + ia y . Therefore (7) holds if and only if a y = 0. 

With respect to an arbitrary fixed basis we can write p{6) = |(1 + u{6) • a) and 
À(^) = 2p(9) = Ü(9) ■ a = r(9)v(9) ■ a where Ü(9) and v(9) are orthogonal unit vectors, 
r{9) = 11^(^)11 is real and nonnegative, and B is the vector of the Pauli spin matrices with 
respect to the fixed basis. Using the familiar relations a 2 x = 1, a x a y = —a y a x = ia z , 
and their cyclic permutations, and the fact that the spin matrices are traceless, one 
finds that the quantum information 1(9) = r(9) 2 . The direction v(9) is uniquely defined 
if r(9) > and from now on we assume this is true for all 9. Then (4) is satisfied 
if and only if for p(x; 9)p(dx) almost all x, m(x) is proportional to a projector for a 
spin direction in the plane spanned by u{9) and u{9). In particular, any simple (von 
Neumann) measurement of spin in a direction in this plane attains equality in (3) if 
both outcomes have positive probability. 

Let C(9) denote the great circle on the unit sphere formed by the intersection of 
the sphere with the plane V(9) spanned by Ü(9) and v(9), let n(9) be the normal unit 
vector to this plane. We suppose that these objects vary smoothly with 9. Recali that 
Ü(9) moves in the direction v(9). Now as 9 varies, either V(9), C(9), and n(9) are all 
fixed or they all vary. In particular, the intersection over all 9 of the planes V(9) is 
either a fixed plane V (equal to V(9) for all 9), or & fixed straight line C, or the origin 
0. In the second case the normal n{9) and the circle C{9) must be rotating about the 
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fixed line C. The rotation must have nonzero speed for vàlues of 9 in a set of positive 
measure. If we assume that the model is identified, so that p(9) is a one-to-one function 
of 9, then at most for two vàlues of 9 can the true spin Ü(9) lie in C. So there exists a 
9 for which the rotation has nonzero speed and Ü(9) does not lie in the axis of rotation. 
But then at this point the derivative of u{9) must have a nonzero component in the 
direction orthogonal to the plane V(9), which is a contradiction. 

So we only have two possiblities: either the plane V{9) is fixed and u(9) is moving on 
the great circle C in the plane, or Ü{9) moves on some other curve and the intersection of 
all V{9) contains only the origin. In the first case any measurement with all components 
proportional to projectors of directions in this plane, and with p(x; 9) > for all x and 
9, achieves the inequality (3) uniformly in 9. Conversely, under the positivity of p(x; 9), 
only such measurements uniformly achieve the bound. In the second case a measurement 
which uniformly achieves the bound would have to have all m(x) equal to zero, which 
is impossible. Thus there is no uniformly attaining measurement in this case. 

For example, consider a spin-half particle in the pure state \ip) = \ip(rj, 9)) given by 

e-^ 2 cos(77/2) \ 

e^sin^/2) )■ { ) 

This pure state has density matrix p == |(1 + u • a) where u = Ü(i], 9) is the point on 
the unit sphere in R 3 with polar coordinates [i], 9). Suppose the colatitude i] e [0,7r] is 
known and exclude the degenerate cases i] = or r\ = n; the longitude 9 G [0, 2n) is the 
unknown parameter. 

We have a pure state so A = 2p = 2Ü ■ a = sin(^) u(n/2, 9 + n/2) ■ a = r(9)v(9). 
The quantum information is r 2 = sin 2 rj. As 9 varies, Ü(9) traces out a great circle if 
and only if r\ = n/2. Consequently, for rj ^ n/2, no measurement M exists with Fisher 
information i(9] M) equal to the quantum information 1(9) whatever the value of the 
unknown parameter 9. If r\ = n/2 it is possible to achieve the bound uniformly in 9. 
Any measurement with everywhere positive density and all components proportional to 
projector matrices for spin directions in the plane r] = n/2 will do the job. A simple 
measurement of spin in one particular direction in that plane attains the information 
bound at all 9 except for 9 equal to that direction or opposite to it. At these points 
the distribution of the outcome is degenerate and the Fisher information not defined. 
However since the Fisher information is continuous (indeed, constant) in 9 this is a 
non-essential singularity. 

4. Asymptotic attainability and vector parameters 

We have shown, for the case of a one-dimensional parameter, that only for rather special 
models will a measurement M exist such that i(9; M) = 1(9) for all parameter vàlues 
9 simultaneously. It is on the other hand possible to find a measurement M such that 
at a given parameter- value, i(9; M) = 1(9), as Braunstein and Caves indicate: take 
each m(x) proportional to a projector onto an eigenspace of the quantum score \(9). 
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They do not remark on the possible dependence of M on 9. However, if all we know is 
that p = p(9) for some 9, we do not know which measurement to use. The eigenspace 
decomposition of À generally depends on 9 so this does not define a measurent M which 
achieves the bound uniformly in 9. This is not the only way to achieve the bound, but 
the previous section shows that one cannot in general expect there to be a uniformly 
attaining measurement. 

Note that the classical information based on n independent and identically 
distributed realisations from a given density p(x, 9) is equal to n times the information 
for one realisation. Similarly, the quantum information in the state p(9)® n corresponding 
to n identical partides each in state p{9) is n times the quantum information for one 
particle. 

Braunstein and Caves' aim was to define a statistical distinguishability mètric 
between quantum states. Suppose the measurement M on a single particle satisfies 
i(9; M) = 1(9). Then the maximum likelihood estimator of 9 based on n separate 
measurements of M on identical copies of the given quantum system, by classical 
results in mathematical statistics, is generally an asymptotically unbiased estimator 
with asymptotic variance (ni(9; M))^ 1 = (nl(9))~ 1 . By the quantum Cramér-Rao 
bound applied to the joint system of n partides, no estimator based on any measurement 
whatsoever on p m (9) can do better. Thus 1(9) appears to exactly characterize the rate 
at which one can determine 9. 

However this argument is flawed since the measurement M involved will be a 
different measurement for each 9, and the whole point is that 9 is not known in advance. 
The question therefore remains: does there exist a measurement procedure not depending 
on 9 on the state p® n , on the basis of which an estimator of 9 can be constructed having 
asymptotic variance (nl(9))~ 1 l If the answer is £ yes', then Braunstein and Caves' 
proposed role for the quantum information 1(9) in defining a statistical distinguishability 
mètric is well motivated. 

It seems rather natural to try a two-stage procedure: first estimate the parameter 
using a perhaps ineflicient procedure on a vanishing proportion of the partides, say 
n = n a (0 < a < 1) out of the total of n; now carry out the 'estimated optimal 
measurement' on the remaining ones. In both stages only simple or von Neumann 
measurements (measurements of classical observables) on separate partides are needed. 

In our example (8) this would reduce to the following. Measure the spin a x on 
k = |n of the copies. The number of +l's observed is binomially distributed with 
parameters k and p = ^(sm rj cos 9 +1). Similarly for another k measurements of the 
spin a y we get a binomial number of ' + 1' with parameters k and p = ^(smi]sm9 + 1). 
This allows us consistent estimation of both siné? and cos 9 and hence of 9 e [0,27r). 
Denote such an estimator by 9. We saw that À in this example was proportional to the 
spin in the direction (n/2, 9 + n/2). Let us use the remaining n' = n — no partides to 
measure this spin with 9 replaced by 9. Given 9, this results in a binomial number X 
of '+1' with parameters n! and p = |(1 — sin?7sm(# — 9)). Let 

9 = 9 + arcsin((n' — 2X)/(n'sin?7)). 
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Analysis of this 'final' estimator shows that 9 has asymptotically the J\f(9, (n sin 2 ^))^ 1 ) 
distribution (the normal distribution with indicated mean and variance), whatever 9, 
so that the quantum ínformation bound is asymptotically achievable by our two stage 
procedure. 

This approach will work in wide generality in problems with a one-dimensional 
parameter 9. Suppose, as typically will be the case, that one can construct a consistent 
estimator 9 based on certain measurements on a vanishing proportion of the partides. 
Compute the quantum score at 9 — 9, and measure it on each of the remaining partides. 
Compute the maximum likelihood estimator 9 of 6 based on the new data, whose 
probability distribution depends on the unknown 9 (as well as on 9, which is at this stage 
fixed). We argue as follows that 9 has approximately the J\í(9, (nl(9))~ 1 ) distribution, 
thus this estimator asymptotically achieves the quantum Ínformation bound. Let i(9; 9) 
denote the Fisher Ínformation for 9 in a measurement, on one particle, of the quantum 
score at 9; thus i(9; 9) = 1(9) for all vàlues of 9, but generally i(9; 9) < 1(9). Now for n 
large, 9 is close to 9. By the classical results for maximum likelihood estimators, given 9, 
9 has approximately the M(9, (ni(9; #)) _1 ) distribution. So if p depends on 9 smoothly 
enough that i(9; 9) is close to i(9; 9) = 1(9) for 9 close to 9, we have that unconditionally 
9 has approximately the M(9, (n/(6>)) _1 ) distribution, hence asymptotically achieves the 
bound. 

Consider now the case of vector parameters. Both quantum and Fisher Ínformation 
numbers are naturally generalized to ínformation matrices; see Helstrom (1976), Holevo 
(1982). The Braunstein and Caves result generalizes to the following result: the 
quantum ínformation matrix is larger (in the sense that the difference is positive semi- 
definite) than the Fisher Ínformation matrix based on the outcome of any measurement 
M. However the bound is no longer attainable. As we saw above, a best measurement for 
each parameter separately is a measurement of the quantum score operator. Typically 
these do not commute and hence cannot be measured simultaneously. On the other 
hand, consideration of all smooth one-dimensional sub-models of a given model p = p(9) 
shows that the quantum Ínformation matrix is the smallest matrix larger than the Fisher 
Ínformation matrix of any measurement on a single particle. 

For instance, suppose we want to simultaneously estimate both parameters i], 9 of 
the pure-state, spin-half system (8); in other words, we have a completely unknown pure 
state. Rename 9 as 0, and let 9 from now denote the vector parameter with elements i], 
<j). Suppose we may dispose of a large number of identical copies of this system. Let 1(9) 
denote the 2x2 quantum Ínformation matrix, and i(9; M) denote the Fisher Ínformation 
matrix based on the outcome of a measurement M, both for a single copy of the quantum 
system. The quantum scores, for a single particle, for the two parameters i] and <f> are 
(J r)+n/2,<f> and sin ?7(T 7r /2,<}í)+7r/2 respectively. After a small proportion of measurements we 
know roughly the location of the parameter, and it is sufncient to investigate optimal 
measurement at a 'known' parameter value. 

Without loss of generality let this be the special point 9 with rj = ir/2, O = 0; 
p — |(1 + a x ). At this point the quantum scores are a y and — a z , and the quantum 
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information matrix is the identity 1. The intersection of the xy and the xz planes 
is the x axis, so it seems that in order to maximize Fisher information both for 77 
and for <fi we should essentially measure spin in the x direction. However the Fisher 
information matrix i{9; M) based on the outcome of this measurement is not defined at 
9 = (n/2, 0) = 6*o since the distribution of the outcome is degenerate. The singularity is 
essential since i{6; M) is not continuous at 6 = 9 . In particular, as one moves towards 
#0 along either of the two great circles formed by varying one of the two components 
of 6, i(9; M) converges to each of the diagonal matrices with diagonal elements 1, and 
0,1. These are the two Fisher information matrices corresponding to von Neumann 
measurements of the two quantum scores, each giving maximal information about the 
corresponding component of 9 and zero information about the other. 

Since different components of the parameter vector have incompatible quantum 
scores, it is clear that for different loss functions, different measurements will be optimal. 
No single procedure will (asymptotically) dominate all others. Moreover, since we 
cannot achieve the quantum information bound by measurements on single partides, it 
is possible that joint measurements on several partides simultaneously could give larger 
Fisher information (per particle) than measurements on separate partides. 

In some very special cases, an optimal procedure is known. An appealing loss 
function in the completely unknown pure spin-half model is one minus the squared 
inner-product between the true state vector and its estimate. This equals one minus the 
squared cosine of half the angle between the points on the Poincaré sphere representing 
the two states. At the special point under consideration therefore, the loss function is 
asymptotically equivalent to one quarter times the sum of the squares of the errors in 
i] and (j). Massar and Popescu (1995), in response to a problem posed by Peres and 
Wootters (1991), exhibited a measurement, optimal in the Bayes sense, with respect to 
this loss function and a uniform prior distribution. It had an asymptotic mean square 
error A/n. This was a genuine generalised measurement of the composite system p® n . 
They showed that for the case of n = 2 there were no measurement methods of the two 
partides separately which were as good as the optimal method, and this is expected to 
hold for all n. 

Instead of this exactly optimal procedure (with respect to the given loss function 
and under a uniform prior) consider taking with probability half measurements of o y 
and a z , independently on each particle. We find that the Fisher information matrix 
(based on one observation) for 77, 0, at 77 = n/2, <p = 0, is |l, or one half of the quantum 
information matrix. The inverse of this matrix, 21/n is an asymptotically achievable 
lower bound to the covariance matrix of (asymptotically unbiased) estimators of 77, <p 
based on n of such measurements. The maximum likelihood method would provide 
estimators asymptotically achieving this bound. The sum of the variances is A/n, the 
same as what is achieved by the Massar and Popescu procedure. 

Thus the following two-stage procedure, similar to what we proposed in the one- 
parameter case, should have asymptotically equivalent covariance matrix to that of the 
Massar and Popescu procedure, and will also be optimal with respect to a uniform 
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prior distribution and any smooth loss function, invariant under rotations of the sphere. 
First carry out measurements of each of a x , a y and o z on a small proportion of separate 
partides. Compute from the results a consistent estimate of (r],(f>). With respect to 
a rotated coordinate system putting the estimated value at (rj, 0) = (n/2,0), measure 
alternately a y and a z on the remaining partides. Estimate (77, 0) in the new coordinate 
system by the method of maximum likelihood using the second stage observations. 
Finally rotate back to the original coordinate system. Note that if we modified this 
scheme by just measuring a x in the second stage, then overall the procedure determines 
the radial distance of 9 from 9 with precision of the order n -1 / 2 , but says nothing about 
the direction, so that finally 9 has not been localized to this precision at all. 

This conjecture has been confirmed by recent further work of Gili and Massar 
(1999). Consider a sequence of measurements on n identical copies of a spin-half state, 
on which is based a sequence of estimators 9. The parameter 9 might be one- or 
two-dimensional for a pure state model, one-, two- or three-dimensional for a mixed 
state. Suppose that the estimators are asymptotically unbiased and have covariance 
matrices asymptotically of the form V(9)/n. Then it is now known that the collection 
of attainable V is precisely {V : trace I(9)~ 1 V(9)~ 1 < 1 for all 9}, in any of the following 
cases: the parameter is one-dimensional, or the state is pure, or the measurement can be 
implemented by separate measurements on separate partides. These optimal limiting 
covariance matrices can all be achieved using a two stage adaptive procedure of the 
type described above. The collection of attainable V corresponds to the collection of 
attainable inverse Fisher information matrices for measurements on single partides. 
Using joint measurements on mixed states with more than one unknown parameter, one 
can attain strictly smaller asymptotic covariance matrices. But a clean description of 
what is attainable is not known. One would like to describe the collections of scaled 
information matrices {i(9; M m ) /m} for m = 1,2,..., where M m is an arbitrary joint 
measurement on m partides. These sets are all convex, they grow with m; we know the 
set for m — 1; and each set is included in the set of matrices less than or equal to 1(9). 
The inverses of these information matrices will be the achievable (scaled) asymptotic 
covariance matrices based on measuring a large number n of partides in groups of m at 
a time. 

To conclude, in the multiparameter case, the bound implied by the quantum 
information matrix is not even asymptotically achievable. The rate at which one can 
distinguish between more than two neighbouring quantum states does not correpond to 
the rate at which one can distinguish between just two; it depends on what aspect of the 
quantum states is important, and it depends on whether one may use joint measurements 
or only separate measurements. The quantum information matrix only plays a role in 
special cases. 
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